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We analyze the developing of bipartite and multipartite entanglement through the Mott-Insulator 
- Superfluid quantum phase transition. Starting from a Mott insulator state, where a filling fac- 
tor v = N/M = 1 per lattice site is considered, we derive an exact expression for a completely 
connected graph configuration of bosons and show how bipartite and multipartite entanglement 
signals the phase transition predicted in previous works. Moreover, through the transition bipartite 
entanglement shows to be monogamous. 
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Entanglement is a central resource of quantum me- 
chanical systems and is particularly important for appli- 
cations in quantum computation and quantum informa- 
tion science [T]. Recently it has been argued that entan- 
glement may be quite relevant in many particle systems 
under quantum phase transitions (QPT) (See eg. [2]). 
Indeed it has been proved that under certain conditions 
a non-analyticity appearing in a many-particle system 
ground state will be signaled in any bipartite [3j and 
multipartite ^ [S] entanglement measure. Several spin- 
1/2 models have been proved to follow this property. On 
the other hand, the interest in bosonic systems has in- 
creased recently due to the actual accessibility to several 
models in optical lattices experiments. One characteris- 
tic model that has attracted considerable attention is the 
Bose-Hubbard (BH) model under a Mott-insulator (MI) 
- Superfluid (SF) QPT. The concepts about phase transi- 
tions in the BH model were firstly derived by Fisher et al. 
[5], and experimentally achieved in a remarkable exper- 
iment by Greiner et al. [5], becoming an intense object 
of investigation [TU] . Some entanglement measures were 
explored previously in a wide variety of physical config- 
urations [HI [131 [m Hi]. At constant densities, an infi- 
nite order Berezinsky-Kosterlitz-Thouless (BKT) quan- 
tum phase transition [17] from the SF to the MI phase 
at low temperatures is expected for this system. One im- 
portant question arises: How is this infinite order QPT 
signaled by entanglement measures? 

In this paper we develop this questioning by investigat- 
ing a typical bosonic system suffering a constant density 
MI- SF QPT. We investigate entanglement in two con- 
figurations through analytical results: (i) A geometry 
referred as completely connected graph (CCG), where 
atomic interactions occur on-site and each atom is al- 
lowed to hop from one to any site of the whole lattice. 
Fig. (l)a. (ii) A linear chain of bosonic atoms trapped 
in a ID optical lattice. Atomic interactions occur on-site 
and each atom is allowed to hop to the two neighboring 
sites. Fig. (l)b. This picture corresponds to the ID BH 
model with mean number of bosons equal to 1. Both 
situations are analyzed considering a fixed filling num- 
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FIG. 1: TV- modes configuration, (a) CCG. (b) BH. 



ber and consequently fixed density, here represented by 
the rate between the number of atoms (N) an sites (M), 
ly = N/M = 1. 

Let us consider the Hamiltonian for bosonic atoms in 
an external trapping potential 

+ ^C/o j dVV'^(r)^t(r)^(r)V-(r), (1) 

where ip{r,t) and tlj^{r,t) represent the field opera- 
tors, V{r) describes the optical lattice potential, Uq = 
m measures the strength of the two-body interac- 
tion, m is the atomic mass and a denotes the s-wave scat- 
tering length. The derivation of the BH or CCG Hamil- 
tonian follows the standard procedures assuming all par- 
ticles to be in the lowest band of the optical lattice [9lfT0]. 
In this treatment, we expand the field operators in terms 
of the Wannicr functions ')jj{x) = '^^aiw{x — Xi), where 
tti is the annihilation operator for a particle in site Xi 
[9l[T^. Then, after some algebraic manipulation, we ob- 
tain the following Hamiltonian using an assumption that 
the total number of atoms A'^ is a conserved quantity 



JY [aja, + a]a,] , (2) 



where U = Airagfi'^ j d'^x\w{x)\'^ [m is the self-collision 
rate or the strength of the on site repulsion of two 
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atoms on the lattice site i and J = J (Pxw*{x — 
Xi) [— (?i^V^) + Vb(a;)] w{x — Xj) is understood as the 
hopping matrix element between adjacent sites z, j or the 
tunneling rate. The operators Ni = a\ai count the num- 
ber of bosonic atoms at the lattice site i; the annihila- 
tion and creation operators and a\ obey the canonical 
commutation relations [ai,a|] — 6ij. The BH Hamil- 
tonian consists of ^ taking the hopping term only be- 
tween neighboring sites, while in the CCG configuration 
we consider hopping between any two sites. Calculations 
through mean field theory, consistent with this last con- 
figuration, indicate the phase transition at the critical 
point oi U/J = 5.6z, where z = 2d is the number of 
nearest neighbors |HJ [TOl |TT]. For a ID optical lat- 
tice, this mean field theory propose a ratio J/U w 0.08 
to the transition point. These results were corroborated 
by the experimental investigation on phase transitions 
in an ID optical lattice developed by Storfele et. al. [7]. 
Moreover there are well-known solutions for its two quan- 
tum phases [11]. Deep into the SF phase the system is 
described by a coherent state where the probability dis- 
tribution for the local occupation of atoms on a single 
lattice site is Poissonian. Furthermore, this state is well 
described by a macroscopic wavefunction with long-range 
phase coherence throughout the lattice 

/ M \^ 

In the MI phase the fluctuations in atom number of a 
Poisson distribution become energetically costly and the 
ground state of the system will instead consist of locahzed 
atomic wavefunctions with a fixed number of atoms per 
lattice site minimizing the interaction energy. The many- 
body ground state is then a product of local Fock states 
for each lattice site [S] 

■ i—l 

Although there is a difference between the two configu- 
rations, BH and CCG, it is known that in the thermody- 
namic limit, TV oo, the behavior of both is very similar. 
Experimentally it has been observed that the systems are 
sufficiently large with a number of sites M of approxi- 
mately 10^ [8]. In this case the SF solution becomes in- 
distinguishable from a coherent state that factorizes into 
a product of local coherent states of energy lattice site. 
This being a consequence of the commutation relation of 
the bosons at different sites [H]. The order parameter 
ATV^ =< A^^ > — < >^ can be used as a signature 
of the transition. In fact, for MI (incompressible) phase 
AA"^ = 0, showing us that the distribution is a Fock one, 
while in the SF (compressible) phase AA^^ 0, and a 
Poissonian distribution is observed. 



The simplest way to investigate the evolution of en- 
tanglement through the QPT is to developed a pertur- 
bation treatment to determine the pure state of the sys- 
tem, departing from the MI ground state Q and taking 
the hopping term as a perturbation [TT]. This so-called 
strong coupling expansion is valid since the transition 
is expected to occur at J/U very small ~ 0.08. Em- 
ploying periodic boundary conditions and considering the 
unperturbed Hamiltonian as Ho = U Ni {N^ — I) and 

yy = ^J2<ij> + '^l^ij where we define A = J/U 

as the perturbation parameter. Since we are dealing 
with the microcanonical ensemble we neglect the chemi- 
cal potential n in our calculations. This is reinforced by 
the assumption that the total density of bosons is fixed. 
We perform our analytical calculations starting from a 
MI state considering an occupation number of exactly 
1/ = N/M = 1 atom per site. We expect thus the MI-SF 
phase transition as soon as A departs from 0. 

Since the system state is pure, multipartite entangle- 
ment (ME) can be detected by the reduced one-site linear 
entropy S = (l ~ '^''Pw i); where d is the dimension 
of the relevant Hilbert space 0. Using the Equations 
([3| and Q is possible to render an analytical expres- 
sion for the reduced (one -site) estate and thus for the 
reduced linear entropy for the CCG configuration, con- 
sidering a general normalized state of the system. For 
the CCG it consists in distributing A^ atoms in M sites, 
considering the occupation number of each site as a label, 
and verify how many states are possible for the system: 
= I'D) = <^i\'i)j where V is the combinatorial form 
to distinctly distribute A^ particles in M ~ N sites, 

X>^(--')^ (5) 

Here a.; is a complex parameter characterizing the contri- 
bution of each combinatorial state \i) in the description of 
the physical state ai is related to the correlations of 
a site with the rest of the lattice. Since in this model we 
deal with only two-body interaction is indeed related 
to two-body correlations, thus bringing all the informa- 
tion about the reduced system. Therefore it certainly is a 
function of the imbalance between the interaction term U 
and the hoping term J. Whether it is continuous or not 
is what determines the type of the phase transition from 
MI to SF. In our approach, a; is obviously a truncated 
power series over the perturbation parameter A. The 
least costing global state starts in the MI phase with one 
single atom occupying each site: |111---1), whose pa- 
rameter is Oi. After some combinatorial reasoning and 
partial trace over A^ — 1 sites the general relation for the 
reduced one-site state for A^ = M is given by 

f{N) N 

p^,i = (A^-l)!^^|a,|^A,,|j)(j|, (6) 

i=l j=0 
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where f{N) is the integer partition function of the num- 



ber N, and Ai 



l/HeU represents all the possibles 



arrangements for the symbols , necessary to write down 
the occupancy j of a particular site. Eq. ^ satisfies the 
constraint determined by \ai\^ > jai+il^, that is under- 
stood as a consequence of the energetic cost since the 
on-site interaction becomes very strong and leads to a 
highly unstable state as the occupation in a determined 
site is increased. Now the reduced one-site linear entropy, 
S, can be calculated from Eq. ^ and plotted in function 
A and N to up to the required order, as in Fig. (2). 

Before we move on to discuss the entropy, let us make 
some remarks. Firstly, we notice that in principle, with 
increasing N it will be required higher orders of pertur- 
bation to describe the exact state. However from the 
hierarchy of above those terms are likely to be very 
small, and thus a first or second order perturbation is 
invariably enough for the description of the behavior of 
ME. Secondly we remark on the possibility of S signaling 
the QPT. Analyzing the first derivative in A of the S we 
find in the thermodynamic limit iV — )■ oo. 



as 

9A 



= -2 



d- 



d\ 



dX 



(7) 



In this limit, one has to consider an infinite occupation 
number and thus an infinite number of a^'s each one of 
the I ail tends to zero and in this same limit Aij is just a 
constant. Thus each i-th term in the summation is close 
to zero unless the derivative in is divergent or at least 
very large. For a continuous QPT it is expected that 
d\ai\'^/dX be the origin of a discontinuity or divergence 
if the transition is of the second order [5^ . For an infinite 
order QPT although the derivative in ai may be large it 
is finite, but still can signal the QPT. The same reasoning 
is valid for derivatives of higher orders, which are then 
continuous. On the other hand the energy corresponding 
to the system perturbed state E{X) = A"£^"\ where 
is the ri-order energy, is always continuous on A as 
well, but it is unlikely to contain any non-analyticity, and 
thus will not signal the QPT. 

In Fig ^ the linear entropy and its derivative are plot- 
ted for a number of atoms up to 1000 for both the CCG 
and the BH configurations, employing a perturbation up 
to second order. To start up, we exemplify in Figs. ([2|a 
and b the two simplest situations, when N = M = 2 and 
3, where the BH and the CCG configuration are indistin- 
guishable. For A small the system can be represented 
as a product of local states, with no entanglement. How- 
ever, as A is increased, we observe an increase of the linear 
entropy, stabilizing after A « 0.3 for iV = 3. This means 
that as the system tends to the SF phase, its state no 
longer is separable, and consequently the entanglement 
is increased. Different behaviors for the BH and CCG 
configurations occur when TV = M = 4 and so on, as can 
be viewed in Figs. (l2|c and d. These distinct evolutions 
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FIG. 2: (a) Reduced one site linear entropy, and (b) corre- 
spondent derivatives, as a function of A for N — A'l = 2 (solid 
line) and 3 (dotted line) for the ID-BH and the GCC config- 
uration, (c) Reduced one site linear entropy, and (d) respec- 
tive derivative, as a function of A for N — M = 4 situation 
in ID-BH (solid line) and CCG (dotted line) configuration. (e) 
Reduced one site linear entropy, and (f) corresponding deriva- 
tive, as a function of A for A'^ = M = 10 (dotted line), 50 
(solid), 100 (triangles) and 1000 (circles) calculated up to the 
first order for the ID-BH configuration. 



are consequence of the different form of the perturbation 
term of the Hamiltonian. In Figs. (|2])e and (|2])f we show 
for the BH model employing first order perturbation only, 
for several distinct number of bosons. The linear entropy 
characterizes very well the two phases predicted by the 
BH model. When A we observe no entanglement, 
characterizing the MI phase. In the other limit, as A ^ 1 
the state tends to be maximally entangled. This is a typ- 
ical behavior of a system in the SF phase. We remark 
that the linear entropy profile is very similar to the or- 
der parameter AA^^ (not shown). Thus AiV^ is indeed 
a witness of the ME. Although we cannot do strong as- 
sumptions about the point where the phase transition 
occurs, given the perturbative approach, the behavior of 
ME signals the QPT correctly. 

We quantify bipartite entanglement (BE) through the 

negativity [19], defined asJ\f{p) = — where 
is the trace norm of the partially transposed state pjj of 
any pair of sites {i, j} for the BH model. This is equiv- 
alent to the absolute value of the sum of negative eigen- 
values of p^^, vanishing for separable states. Although 
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FIG. 3: (a) N12 (dashed line) for iV = M = 3 and the reduced 
one site hnear entropy (sohd hne) and (b)A/i2 (dashed Une), 
A/13 (dotted hne) and the reduced one site hnear entropy for 
N = M = 4. 



strictly necessary and sufficient only for Hilbert space di- 
mension up to 2 (g) 3, in the present case it is correctly 
characterizing the BE for any A^. In Fig. ([3| we plot the 
negativity and the linear entropy for = 3 and = 4. 
Since for iV = 3 there are only nearest neighbor for any 
site the negativity is the same for any two sites. For 
= M = 4 the negativity shows different behavior for 
nearest and next-nearest neighbors, here exemplified by 
P12 and P13, respectively. Contrary to spin-1/2 model 
(see e.g. [12]) the next-nearest neighbor negativity can- 
not be neglected in comparison to the nearest neighbors 
one. Instead they show a peculiar behavior signaling the 
monogamy of entanglement. In the simplest situation 
(A^ = M = 2), since the system is pure, the negativity 
is equal to the one site reduced linear entropy as there is 
only BE. When N = M = 3 the negativity after showing 
an increase in the transition, stabilizes at lower values 
than the one site linear entropy. This is a signature that 
genuine tripartite entanglement do exist in the system. 
Indeed the ME develops through the bipartite one (the 
two curves evolve similarly for small A). But as soon as 
tripartite entanglement starts to develop the two curves 
diverge. In the N = M = 4 case, the negativity for pi2 
and yOi3 shows that distinct kind of BE exists. A/12, after 
increasing with the entropy, decreases as A/13 increases 
(see the inset in Fig. 2), showing an interesting bound 
on the distribution of BE. Subsequently both measures 
stabilize to closer values. This bound on BE is typical of 
the monogamy of entanglement [T5] in the sense that the 
increasing of the pi2 entanglement at the same interval 
where the yOi3 one decreases and a reciprocity maintain a 
constant amount of the BE. Here the difference between 
the ME and the two types of BE, namely the residual 
entanglement signals both tripartite and quadripartite 
entanglement. This feature continues for increasing A'^, 
and there is an increasing number of types of BE. The 
initial crossing of BE decreases with N, stabilizing after a 
while at higher values. This is an evident manifestation 
of the many classes of entanglement present in the BH 
model. The SF phase is intermediated by a strong ME 
state where any mode (site) is entangled with the others 



in many distinct ways. 

In short, we analyzed ME and BE for the CCG and 
ID-BH configuration of bosonic atoms trapped in an 
optical lattice. We described the behavior of entangle- 
ment through the MI-SF QPT predicted for these models 
showing how it is signalized by two entanglement mea- 
sures, the reduced one site linear entropy and negativity. 
By employing the linear entropy just the diagonal ele- 
ments of the whole density matrix must be taken into 
account, since the partial trace keeps only these terms 
(see Eq. (|6|). It represents a remarkable reduction of 
numerical resources for investigation of a QPT in com- 
parison with the full Hamiltonian diagonalization. For 
instance, for A^ = M = 10, by Eq. (|5| one must diago- 
nalize a 92378 x 92378 matrix, while through the linear 
entropy it involves the calculation of just 42 coefficients 
in the case of GCC configuration, which can be solved by 
some computational method with relative simplicity. 
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